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Abstract 

We derive and then enhance a generalized matrix version of Pellet's theorem, itself 
based on a generalized Rouche theorem for matrix- valued functions to generate upper, 
'• lower, and internal bounds on the eigenvalues of matrix polynomials. The matrix 

C/^ [ version of Pellet's theorem is subsequently used to improve the original scalar theorem. 
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1 Introduction 

Polynomial eigenvalue problems have been investigated for quite some time ([7], [12], |13j ) 
and have important applications in a wide range of engineering fields such as vibration 
analysis, acoustics, and fluid mechanics - to name just a few (|19j). It is, in general, costly 
to compute polynomial eigenvalues for large problems, but bounds on such eigenvalues are 
relatively easy to obtain. They provide information on the location of eigenvalues that can 
be used by iterative methods for computing them and are also useful in the computation 
of pseudospectra. 

The polynomial eigenvalue problem is to find a nonzero eigenvector v, corresponding 
to an eigenvalue z satisfying P{z)v = 0, where 

P{Z) = AnZ^ + + • • • + Ao, 

with Aj £ ^rnxm j _ . . . n. We will assume throughout that det (P(z)) is not 
identically zero. If An is singular then P has infinite eigenvalues and if Aq is singular then 
zero is an eigenvalue. There are nm eigenvalues, including possibly infinite ones. The 
finite eigenvalues are the solutions of det {P{z)) = 0. 

To explain the aims of this work, we start with the scalar polynomial with complex 
coefficients 

p{z) = anz"' + a„_iz""^ H h ao , 

and a^aQ ^ 0. If C is a zero of p, then for ^ with l<A;<n— Iwe have 

-akC'' = anC + an-iC"^^ + • • • + flfc+iC''^^ + a^-iC^'^^ + ■ ■ ■ + uq , 
and an analogous equality for k = 0,n, which implies that 

lofcllCl'' = \anC + an-iC~^ + ■ ■ ■ + ak+iC''^^ + Uk-iC'"'^ + ■ ■ ■ + ao 

< |a„| iCr + l«n-i| ICr' + ■■■ + \ak+i\ |Cl'+' + \ak-i\ ICI'-' + • ■ • + |ao| •(!) 
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Applying the equivalent of inequality ([T]) with k = n, we obtain that \C\ must satisfy 



\an\\Cr - K-iWCr' |ai||C|-|ao| <0. 

This means that \(\ can be no larger than the unique positive root of 

lonlx*^ — \an-i\x"'~^ — . . . — \ai\x — \ao\ = . 
Analogously, one finds that |C| can be no smaller than the unique positive root of 

jonlx"" + \an-i\x^~^ + • • • + |ai|a; — |ao| = . 

These results are due to Cauchy ([3], [HI Theorem (27,1), p. 122]). 

When we apply inequality ([1]) with 1 < k < n — 1 and 7^ 0, we obtain 

\an\ ICr + K-i\ iCr-' + • • • + l«fc+i| |C|'+' - \ak\ Id' + la^-il ICI'"' + • • • + |ao| > . (2) 
Now, the equation 

la^lx"- + lon-ilx"-"^ H h \ak+i\x'''^^ - \ak\x'' + |afc_i|j;^"^ H h |ao| = 

has, by Descartes' rule, either two or no positive roots. If it has no positive roots or if 
the two positive roots coincide, then inequality ([2]) is always satisfied and provides no 
information on |d- However, if it has two distinct positive roots, say xi and X2 with 
xi < X2, then inequality ([2]) implies that either |C| < xi or \(\ > X2 and that no zero of p 
has a modulus in (xi, 3:2). As a straightforward consequence of Rouche's theorem, Pellet's 
theorem makes this result more precise. It is stated as follows. 

Theorem 1.1 (T'ellet, [I^/, [I^ Th.(28,l), p. 128]) Given the polynomial p{z) = ^ 
an-iz'^^^ + - ■ ■+aiz+aQ with complex coefficients, 0^0^ 7^ 0, andn > 3. Let 1 <k < n — 1, 
and let the polynomial 

fk{x) = + |a„_i|a;''-^ + • • • + |afc+i|x''+^ - |afc|x'' + |afc_i|a;^'^^ + • • • + |ao| 

have two distinct positive roots xi and X2 with xi < X2. Then p has exactly k zeros in or 
on the circle \z\ = xi and no zeros in the annular ring xi < \z\ < X2- 

This means that Cauchy's results can be considered as a special limit case of Pellet's 
theorem. 

We first generalize Pellet's theorem (and therefore also Cauchy's results) to matrix 
polynomials. The standard proof of Pellet's theorem uses Rouche's theorem for analytical 
functions, which suggests that a generalized version of Rouche's theorem for analytical 
matrix-valued functions will be needed. Such a theorem can be derived as a special case of 
the generalized Rouche theorem for bounded linear operators in [HJ Theorem 9.2, p. 206] 
(see also [6]). It is then a fairly straightforward exercise to generalize Pellet's theorem. 
We note that, based on a different and slightly more restricted generalization of Rouche's 
theorem, a generalized Pellet theorem was recently also derived in [3], where it was used 
to determine initial approximations for the numerical computation of the eigenvalues of 
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a matrix polynomial with iterative methods like the Ehrlich-Aberth method ([T], [2], [S]). 
To our knowledge, [3] is the first time Pellet's theorem has been generalized to matrix 
polynomials. 

Our main goal is to find ways to overcome situations where in Theorem 11.11 (or an 
analogous polynomial for the generalized version of that theorem) does not have positive 
roots for a particular value of k. This problem does not seem to have been addressed else- 
where in the literature. It occurs frequently even when there is a significant gap between 
groups of zeros (or eigenvalues). In addition, only polynomials with very special coeffi- 
cients would have more than a few values of k for which the theorem can be applied. We 
will explore two ideas to tackle this problem: one for both scalar and matrix polynomials 
and one for scalar polynomials, which in some cases will lead, perhaps paradoxically, to 
an improvement of the original (scalar) version of Pellet's theorem by using its generalized 
(matrix) version. At the same time, this approach can also improve the upper and lower 
bounds resulting from Cauchy's result mentioned above. 

We will not dwell on the numerical problem of finding the roots of functions like fk in 
Theorem 11.11 if they exist. An efficient way to compute them can be found in jl5j . while 
bounds on such roots were derived in [3]. The number of arithmetic operations that this 
requires is typically dwarfed by the much more costly computation of the eigenvalues of a 
matrix polynomial. Moreover, this work focuses on the enhancement of Pellet's theorem 
if and when this theorem is being used. It does not focus on whether or not it should be 
used in the first place, as this depends very much on external factors. 

The organization of the paper is as follows. In Section 2 we review a few preliminaries 
that will be needed in the later sections. In Section 3 we derive the generalized Pellet 
theorem and in Section 4 we propose two strategies to improve it, along with numerical 
examples. 

2 Preliminaries 

Since it is mentioned several times, we begin by stating Rouche's well-known theorem. 

Theorem 2.1 (Tlouche, !118^ . I 111 Theorem 1-6]) Let f and g be analytic in the interior 
of a simple closed curve C and continuous on C , and let \g{z)\ < \f{z)\ for all z & C. 
Then f + g and f have the same number of zeros in the interior of the curve. 

All matrix norms throughout this work are assumed to be vector-induced (or subordinate). 
The main matrix norms we will use are the 1-Norm, oo-Norm, and the 2-Norm (or spectral 
norm), defined p. 294-295]) for A G (T"''" with elements a^j by 
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The zeros of a monic polynomial p{z j = z" + an-iz 
are the eigenvalues of its n x n companion matrix 



' + with complex coefficients 
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Likewise, the eigenvalues of the monic matrix polynomial 

P{z) = Iz"^ + An-iz""-^ + ---+AiZ + Aq 



where Aj g (T'^x"^ 
companion matrix 



for j = 0, . . . , n — 1, are given by the eigenvalues of the nm x nm block 



C{P) 



-^0 \ 
-Ai 



\ I -An-lJ 

with I the mxm identity matrix. Since the size of / will usually be clear from the context, 
it will be omitted from the notation. 

The reciprocal polynomial of p{z) = + a„„iz"~^ + • • • + ao with oq 7^ is defined 
by Pr{z) = z^p{l/ z) / aQ. Its zeros are the reciprocals of the zeros of p. Likewise, the 
reciprocal matrix polynomial of P(z) = I z^ + An-iz"^'^ + • • • + Aiz + Aq with. Aj S C™^"^ 
for J = 0, . . . , n — 1 and ^0 nonsingular is defined by Pr{z) = z'^Aq^P{1/z). Its eigenvalues 
are the reciprocals of the eigenvalues of P. 



3 Generalized Rouche and Pellet theorems 

To formulate the operator generalization of Rouche's theorem that we will use later, we 
need the following definitions for which we have used similar notation and style as in (8l 
Chapter XI]. We denote by C{X,Y) the space of all bounded linear operators from X 
to Y, where X and Y are complex Banach spaces. When X = Y, this space is denoted 
by C{X). A bounded linear operator S : X ^ Y, where X and Y are complex Banach 
spaces is called a Fredholm operator if its range Im S is closed and if dim(Ker S) and 
dim(y/ImS') are finite. If X and Y are finite dimensional, then S is always Fredholm. 

For an open connected subset of (T, let Q : J] — C{X) be a bounded operator 
function that is analytic on Vt and let F be a simple closed curve in Vt such that its inner 
domain is a subset of Vt. Then Q is said to be normal with respect to F if Q{z) is invertible 
for all z S F and Q{z) is Fredholm on the inner domain of F. Finally, we recall that a 
Hilbert space is separable if and only if it has a countable orthonormal basis. 

The following theorem generalizes Rouche's theorem for operator-valued functions. We 
have stated it as it appears in jH p. 206], but a similar version was already published in 
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Theorem 3.1 (Generalized Rouche theorem for operators.) (J^, J^) Let W,S : 
— C{H) be analytic operator functions, where S7 is an open connected subset of(C and 
H is a separable Hilbert space, and assume that W is normal with respect to the simple 
closed curve F C $7. 

// < 1 for all z £ T, then W + S is also normal with respect to T and 

W + S and W have the same number of finite eigenvalues inside T, counting multiplicities. 

The norm used in the previous theorem can be any norm, induced by a norm on H. We 
obtain the following generalized Rouche theorem for analytic matrix-valued functions as 
an immediate consequence of Theorem 13.11 We note that matrix polynomials are special 
cases of such functions. 

Theorem 3.2 (Generalized Rouche theorem for matrices.) Let A,B:^} _^([;mxm 
be analytic matrix-valued functions, where 0, is an open connected subset ofdJ and assume 
that A{z) is nonsingular for all z on the simple closed curve F C 17. 

// 1 |yl(2;)~"'^i?(2;)| I < 1 for all z S V, t/ien det(A + i?) anc? det(yl) have the same number 
of zeros inside T, counting multiplicities. 

Proof. The proof follows directly from Theorem 13.11 In this case, H = (T", a finite 
dimensional Hilbert space which is (trivially) separable and C (H) = Further- 
more, since A is a matrix, it is a bounded linear operator that is normal with respect 
to F because A(z) is invertible for all z S F and, since (T" is finite dimensional, A{z) is 
(trivially) a Fredholm operator for any z in the interior of F. By Theorem 13.11 A + B is 
then also normal with respect to F and has the same number of finite eigenvalues inside 
F as A, which means that det(A + B) and det(74) have the same number of zeros inside 
F, counting multiplicities. O 

Remarks. (1) The norm in Theorem 13. 21 can be any matrix norm, induced by a norm 
on (T™. 

(2) Since \\A"\z)B{z)\\ < \\A^\z)\\ ■ \\B{z)\\, it is sufficient that \\B{z)\\ < \\A~^{z)\\-^ 
for the condition ||74~-'^(2;)i3(z)|| < 1 to be satisfied. 

Finally, Theorem 13.21 leads to the following generalization of Pellet's theorem (Theo- 
rem [TTT|). 

Theorem 3.3 (Generalized Pellet theorem.) Let 

P{z) = AnZ^ + An-lZ^~^ + ---+AiZ + Aq 

be a matrix polynomial with n > 2, Aj G (U'^XTn j _ _ _ _ ^ Q^^f]^ ^ q_ j^^^ 
invertible for some k with 1 < k < n — 1, and let the polynomial 

fk{x) = \\An\\x^+\\An-i\\x^-'+- ■ ■+Uk+i\\x''+'-\\A~^\\-'x''+\\Ak-i\\x''-'+- ■ •+||^i||x+||^o 

have two distinct positive roots xi and X2 with xi < X2. Then P has exactly km zeros in 
or on the disk \z\ = xi and no zeros in the annular ring xi < \z\ < X2. 
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Proof. By Theorem if 



Akz'' ] ^ f A„z" + An-iz"" ^ + • • • + Ak+iz''^^ + Ak-iz^ ^ + ■ ■ ■ + Aiz + Aq 



< 1 



(3) 

for \z\ = X, then P and A^z^ both have the same number of eigenvalues in the open disk 
1^1 < X, namely, km. 

From Remark (2) following the proof of Theorem [321 ^ sufficient condition for inequal- 
ity ([3]) to be satisfied is 



A„z" + An-iz"" ^ + • • • + Ak+iz^^^ + Ak-iz" ^ + --- + Aiz + Ao 
Since 



< 



AkZ 



(4) 



A„z" + An-iz""-^ + • • • + Ak+iz''+^ + Ak^iz''-^ + ■ ■ ■ + Alz + Ao 
< \\An\\ l^l" + \\An-i\\ 1^1""^ + • • • + ll^fc+ill + \\Ak-i\\ \z\''-^ + ■■■ + \\Ai\\ \z\ + 

inequality will be satisfied for = x whenever x is such that 

(5) 

i.e., when fk{x) < 0, where was defined in the statement of the theorem. By Descartes' 
rule, fk has either two or no positive roots. Since it was assumed that fk has two positive 
roots xi and X2 and since /fc(0) > 0, inequality ([5]) will be satisfied for any x such that 
xi < X < X2. This concludes the proof. O 

Remarks. (1) Applying Theorem lS.SI to the matrix polynomial A'f^^P yields better values 
for xi and X2 since ||A^^Aj|| < ||A^^|| ||Aj||. The disadvantage is that all the coefficients 
need to be multiplied by A'^^ , which could be costly. 

(2) A version of Theorem 13.31 (applied to A'^^P) was also obtained by Bini, Noferini, 
and Sharify in [3] by using a different generalization of Rouche's theorem from |16j and 
|20| . which limits their result to the spectral norm. Because it is based on Theorem 13.21 
Theorem 13.31 holds for any matrix norm induced by a vector norm. This becomes more 
important as the size of the coefficient matrices increases because the spectral norm can 
be costly to evaluate, whereas, e.g., the 1-Norm and oo-Norm can be computed cheaply. 

(3) Whether the polynomial fk will have two positive roots depends on the relative sizes 
of its coefficients. Generally speaking, this will happen when the coefficient of x^ is 
large compared to most of the other coefficients, but it would be very difficult to predict 
how large it needs to be. As in the scalar case, if fk does not have positive zeros, then 
Theorem 13.31 provides no information about a possible gap between groups of eigenvalues. 
In the next section we propose ways to sometimes overcome this problem. 

The following theorem, which is a generalization of Cauchy's result, can be considered 
as a special limit case of the generalized Pellet theorem just as in the scalar case. Its proof 
is entirely analogous and will therefore be omitted. 



Theorem 3.4 (Generalized Cauchy theorem.) All eigenvalues of the matrix polyno- 
mial 

P{z) = Anz^ + A^_^z^-^ + ... + Aiz + Ao, 

where Aj G (urnxm ^ j _ q _ ^ ^^^ /^g j^;] < n when An is nonsingular, and lie in 
\z\>r when Aq is nonsingular, where R and r are the unique positive roots of 

IIA-^ir^x" - WAn-iWx""-^ \\Ai\\x - \\Ao\\ = 

and 

WAnWx"" + WAn-lWx""-^ + ■■■ + \\Ai\\x - 11^0 ^ir^ = ' 

respectively. 

This theorem was also derived in a different way as Lemma 3.1 in [9]. Again, it can 
be apphed to A~^P and A^ P to yield tighter bounds, but at the expense of adding n 
multiplications of two m x m matrices. 



4 Enhancements of the generahzed Cauchy and Pellet the- 
orems 

4.1 Squaring and repartitioning 

The following lemma shows that by squaring the block companion matrix and repartition- 
ing the result, additional Cauchy-type and Pellet-type bounds can be generated. 

Lemma 4.1 The squares of the eigenvalues of the monic matrix polynomial 



P{z) = Iz^ + A„_iz"-^ + --- + Aiz + Aq, 

where n is a positive even integer and Aj G (27™^"^ for j = 0, . . . , n — 1, are given by the 
eigenvalues of the matrix polynomial 



Q{z) = Iz'^/'^ + S 



n/2-l 



.n/2-l 



+ --- + Biz + Bo 



where Bj G (T 



2mx2m 



for J = 0, . . . , n/2 — 1, with 



A. 



-Au^iAq 



Ax -An^iAi + A 



; Bj 



A 



2i 



A 



2j+l 



-An-lA2j + ^2j-l 
-An-lA2j+l + A2j 



i = l,2,...,n/2-l . 



Proof. The companion matrix of P is defined by 



C{P) 



/o 



-^0 \ 
-Ai 



-An-l) 
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so that C'^{P), whose eigenvalues are the squares of those of P, is given by 



/o 



/ 

V 



-Ao An-iAo \ 
-Ai An-iAi - Aq 

-A2 An-lA2 - Ai 



I -A 



n-l 



42 



An-2) 



Since n is even, C'^{P) can be repartitioned into n/2 blocks of size 2m x 2m as follows: 
















/ 








/ 



-^0 

-Ai 


An-lAo 
An^lAi - Ao 




-A2 
-As 


An-lA2 - Ai 
An-lAs - A2 



\ 





I 










/ 



-An-2 An-\An-2 — A^-Z 
-An— I ^n— 1 ~ An— 2 



In other words, C'^{P) can be written as 



C72(P) 



/o 



where 



^0 
Ai 



-An-lAo 
-An-lAi + Aq 



A 



A 



2j 



2j+l 



-Bo \ 
-Bi 



-Bn/2-1/ 

-An-lA2j + A2j-l 
-An-lA2j+l + A2j 



1,2, 



,n/2-l , 



and where / now stands for the 2m x 2m identity matrix. The expression for C'^{P) in 
is the block companion matrix C{Q) of the matrix polynomial 



Q{z) 



,n/2 



+ Bn/2-lZ + 



+ Biz + Bo 



(6) 



whose eigenvalues, being the same as those of C'^{P), are the squares of the eigenvalues 
of P. D 

We can now apply Theorem 13.31 and Theorem 13.41 to obtain information about the 
eigenvalues of P by using C{Q) with Q as in Q that can potentially improve results 
obtained for C{P). This may take the form of improved upper and/or lower bounds, but 
it may also be the case that the conditions for the application of Theorem 13.31 can now be 
satisfied, whereas before they could not, so that a gap between groups of eigenvalues can 
be computed that would otherwise not have been detected. Furthermore, the degree of Q 
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is only half that of P which has a generally beneficial effect on computations involving it. 
Of course, the cost of squaring C{P), while still orders of magnitude lower than that of 
actually computing the eigenvalues, may not be negligable, depending on the size of m. 

Additional bounds can be derived by applying Lemma 14.11 to the reciprocal matrix 
polynomial Pr, which yields a matrix polynomial that we designate by Qji. This matrix 
polynomial is, in general, different from Qr which is the reciprocal polynomial of Q, defined 
in ([6]). The reciprocal of the Cauchy upper bound for P^ or Qr is equal to the Cauchy lower 
bound for P or Q, respectively. However, the reciprocal of the Cauchy upper bound for Q/j, 
while being a lower bound on the moduli of the eigenvalues of P, is, in general, different 
from the Cauchy lower bound for Q. A analogous situation exists for the reciprocal of the 
Cauchy lower bound for Qr. 

When n is odd we consider zP{z) instead of P{z), which simply has m extra zero 
eigenvalues. In such a case, is the zero matrix, which makes Bq in Lemma [4.11 singular 
so that Q cannot be used to obtain a lower bound. To remedy this, one can use zPr{z): the 
reciprocal of the upper bound on its largest eigenvalue provides the desired lower bound 
on the smallest eigenvalue of P. 

Many different additional bounds can be derived by considering MP or NQ instead 
of P or Q for any nonsingular matrices M and N, or by using different scalings of the 
companion matrices. It would therefore be impractical to compare all possible variants. It 
is also difhcult if not impossible to predict which particular bound will outperform other 
bounds for any given situation. Instead, we will present several numerical examples to 
illustrate that our bounds can be valuable complements to existing ones. 

Our bounds' complexity for a particular value of k is linear in the degree n and further 
depends on the coefficient matrices' properties. If, for instance, the coefficients exhibit 
structure, such as symmetry or sparsity, then the various matrix manipulations involved 
in their computation generally require fewer operations. The choice of bounds therefore 
ultimately depends on the situation, which is outside the scope of our discussion. 

We observe that Lemma |4. II can also be applied to scalar polynomials, leading to the 
slightly counterintuitive situation where the generalized matrix version of the theorem is 
used to improve the original scalar version. 

Lemma 14.11 always provides upper and lower bounds when used together with The- 
orem 13. 4[ but it does have limitations when used with Theorem 13.31 since the 2m x 2m 
matrix coefficients of Q are twice the size of those of P. This leads to the computation of 
gaps between groups of eigenvalues that contain an even multiple of m, instead of just a 
multiple of m. Although it is easy to find matrix polynomials for which Lemma |4. II gives 
worse bounds, it is just as easy to find examples where it delivers better bounds, as will 
be shown in the examples below. 

Finally, we remark that higher powers of C{P) could similarly be used, although the 
computational cost involved may render the resulting bounds inefficient. 

In Example 1 we compare the upper and lower bounds resulting from Theorem 13.41 
combined with Lemma 14.11 for both scalar and matrix polynomials. In Example 2 and 
Example 3 we compare the bounds resulting from Theorem 13.31 combined with Lemma [4. II 
for matrix polynomials and in Example 4 and Example 5 we do the same for scalar 
polynomials. These examples focus somewhat more on Theorem 13.31 than on Theorem 13.41 
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since the kind of bounds resulting from it are less frequently mentioned in the literature. 
Although n is an even number in the examples, this is not essential and similar conclusions 
can be drawn for odd numbers. 

Example 1 



For this example we have generated matrix polynomials of degree n = 10 with m = 
1, 2, 10, 25, where m = 1 corresponds to a scalar polynomial. They take the form 

P{z) = Iz^^ + AgZ^ + --- + Ao, 

where each of the m x m matrices Aj, for j = 0, . . . , 9, has complex elements whose real 
and imaginary parts are uniformly distributed in [—1, 1], multiplied by a random number 
which is uniformly distributed in [0, 10]. To each P corresponds a polynomial Q, given by 

Q{z) = Iz^ + B^z^ + --- + Bo, 

where the 2m x 2m matrices Bj, for j = 0, ...,4, are defined as in Lemma 14. li For 
each value of m, one thousand such random matrix polynomials were generated. We then 
compared the ratios of the bounds to the moduli of their largest and smallest eigenvalues, 
respectively. The following upper and lower bounds were compared: 

Upper bounds: P : Theorem 13.41 applied to P 

Q : Theorem 13.41 applied to Q with Lemma 14.11 

Lower bounds: P : Lower bound in Theorem 13.41 applied to P 

Q : Lower bound in Theorem 13.41 applied to Q with Lemma 14.11 
Aq^P : Lower bound in Theorem 13.41 applied to Aq^P 
Bq^Q : Lower bound in Theorem 13.41 applied to Bq^Q with Lemma l4.ll 
Qr : Reciprocal of upper bound in Theorem 13.41 applied to Qr. 

Preceding each bound's description is its corresponding designation in the tables below. 
We recall that the bound designated by Qr is obtained by forming Qr from Pr just as Q 
is formed from P. 

In Tabled] we have listed as percentages the mean of the ratios of the upper bound to 
the modulus of the largest eigenvalue with their standard deviation between parentheses. 
For instance, a ratio of 3/2 corresponds to 150%. The closer the number is to 100, the 
better it is, both for upper and lower bounds. The bounds were computed with the 1- 
Norm, oo-Norm, and the 2-Norm. In Table [2] we have listed the number of times each 
bound was the better bound. In Table [3] and Table H] the same was done for the lower 
bounds. 

As Table [T] shows, using Lemma 14.11 produces better upper bounds on average for all 
three norms, except for the case of scalar polynomials. This advantage seems to increase 
with increasing m. The standard deviations are generally comparable, except for the 
higher values of m with the 2-Norm, where they are larger for the Q bound. The 2-Norm 
gives the best results but is, in general, costly to compute. It is also clear from Table [2] 
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that even when a bound is worse on average, there is still a non-negligible number of cases 
where that bound prevails over the other one. 

For the lower bounds, Table [3] and Table H] show that the advantage goes to for 
m > 1, while P dominates for m = 1. As expected, the bounds Aq^P and Bq^Q are better 
than P and Q, respectively. There were no instances where the latter two delivered the 
best bound. As for the upper bounds, the 2-Norm gives the best results. The advantage 
of using Lemma l4. II is. as with the upper bounds, more pronounced for higher values of 
m. 



m 


P 


Q 


m 


P 


Q 


m 


P Q 


1 


116 (17) 


131 (18) 


1 


116 (17) 


122 (18) 


1 


116 (17) 


118 (16) 


2 


161 (32) 


162 (29) 


2 


162 (33) 


147 (30) 


2 


140 (28) 


133 (23) 


10 


316 (45) 


247 (47) 


10 


317 (46) 


216 (50) 


10 


167 (19) 


161 (31) 


25 


481 (62) 


312 (67) 


25 


481 (62) 


271 (71) 


25 


177 (14) 


174 (39) 



1-Norm 



oo-Norm 



2-Norm 



Table 1: Comparison of the upper bounds with n = 10. 



m 


P 


Q 


m 


P 


Q 


m 


P 


Q 


1 


939 


61 


1 


759 


241 


1 


557 


443 


2 


583 


417 


2 


275 


725 


2 


342 


658 


10 


152 


848 


10 


136 


864 


10 


317 


683 


25 


93 


907 


25 


90 


910 


25 


340 


660 



1-Norm 



oo-Norm 



2-Norm 



Table 2: Best bound frequencies for the upper bounds with n=10. 



m 


P 




Q 






Qr 


1 


87 (9) 


33 


(21) 


87 (9) 


43 (24) 


81 (9) 


2 


49 (13) 


16 


(13) 


65 (13) 


25 (18) 


67 (11) 


10 


10 (6) 


4 


(4) 


27 (10) 


8 (10) 


39 (9) 


25 


4(3) 


2 


(2) 


16 (7) 


5 (6) 


29 (8) 



1-Norm 



m 


P 


Q 


A^'P 




Qr 


1 


87 (9) 


34 (21) 


87 (9) 


43 (24) 


83 (9) 


2 


49 (13) 


16 (13) 


65 (13) 


26 (18) 


70 (11) 


10 


9(6) 


4(5) 


27 (10) 


9 (10) 


42 (10) 


25 


4(3) 


2(3) 


16 (7) 


5(6) 


31 (9) 






oo- 


Norm 







m 


P 


Q 


A-'P 




Qr 


1 


87 (9) 


38 (23) 


87 (9) 


45 (25) 


87 (8) 


2 


61 (14) 


21 (17) 


71 (13) 


29 (20) 


78 (12) 


10 


25 (10) 


7(9) 


39 (12) 


11 (13) 


54 (12) 


25 


15 (7) 


4(6) 


26 (9) 


6(9) 


42 (11) 








2-Norm 







Table 3: Comparison of the lower bounds with n = 10. 
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m 






Qb. 


m 






Qr 


m 






Qr 


1 


761 


9 


230 


1 


720 


22 


258 


1 


394 


15 


591 


2 


395 


19 


586 


2 


277 


19 


704 


2 


175 


11 


814 


10 


85 


17 


898 


10 


99 


7 


894 


10 


76 


9 


915 


25 


48 


11 


941 


25 


56 


3 


941 


25 


61 


3 


936 



1-Norm 



oo-Norm 



2-Norm 



Table 4: Best bound frequencies for the lower bounds with n = 10. 



Example 2 

In this example we illustrate the effect of Lemma 14.11 on the ability of Theorem 13.31 to 
produce a gap between the moduli of two groups of eigenvalues. To do so we have gener- 
ated two differents classes of matrix polynomials of degree n = 14 with ?ti x m = 25 x 25 
matrix coefficients of the form 

P{z) = Iz^^ + yli3zi3 + ... + Ao. 

The elements of the matrices Aj in the first class have real and imaginary parts that are 
uniformly distributed in [-75, 75], [-250, 250], [-5, 5], and [-1, 1] for j = 7,j = 8, j = 9, 
and j 7^ 7, 8, 9, 13, respectively. One thousand such polynomials were generated for which 
the real and imaginary parts of the elements of A13 are in the intervals [—77,17], where 
77 = 0.25,0.125,0. 

For the second class, the elements of the matrices Aj have real and imaginary parts 
that are uniformly distributed in [—75,75], [—250,250], and [—1,1] for j = 11, j = 12, 
and j 7^ 11, 12, 13, respectively. One thousand such polynomials were generated for which 
the real and imaginary parts of the elements of A13 are in the intervals [—77,77], where 
r] = 0.5,0.25,0.125,0. 

For each matrix polynomial, we checked if Theorem 13.31 was applicable with and with- 
out the use of Lemma [4.1l and if it was, computed the ratio of the gap between the moduli 
of the two groups of eigenvalues derived from that theorem to the actual gap. The theorem 
was applied to both P and Q, and also to A^^P and BJ^^^Q, with k = 8 for the first class 
and k = 12 for the second class of matrix polynomials. 

Throughout this example, only the 1-Norm was used. Similar results are obtained for 
the oo-Norm. Although it generally delivers better results, we considered the 2-Norm's 
computational cost to be too high to make it practical. 

On the left in Table [5] are listed, for the first class of polynomials, the number of times 
(out of 1000 cases) Theorem 13 . 31 produced a gap between the mk = 200 smallest and m{n — 
k) = 150 largest eigenvalues for each value of r/ when applied to the matrix polynomials 
A'^^P and B^^^Q with k = 8. The table contains the total number of computed gaps 
for each matrix polynomial as well as the number of times each matrix polynomial was 
the only one of the two for which a gap could be computed. On the right in Table [5] are 
listed, as percentages, the means of the ratios of the gap from Theorem 13.31 to the actual 
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gap, with the standard deviations between parentheses. The rightmost column hsts the 
percentage of cases where the gap for BJ^^^Q was larger than the one for A'j^^P when a gap 
was produced for both. Not enough gaps were computed when Theorem 13.31 was applied 
to P and Q to make a comparison worthwile. 

Tableland Table [7] are the analogous tables for the second class of matrix polynomials 
with k = 12. Here, we have also compared the results obtained by applying Theorem l3.3l to 
P and Q. Overall, the results obtained by applying Lemma [4.11 improve as H^isH becomes 
smaller, which is understandable since it tends to keep \\Bj\\ (j = 0, . . . ,n/2 — 1) of the 
same order of magnitude as (j = 0, . . . , n — 1). This is more important here than 

for mere upper and lower bounds as in Example 1, because the existence or nonexistence 
of positive roots for fk in Theorem 13.31 is heavily dependent on the relative magnitudes of 
its coefficients. 

The results for A; = 8 in Table [5] show clearly that the application of Lemma 14.11 
improves the number of times Theorem 13.31 can be applied. Although it does not always 
increase the average gap width, it was observed that the gap for B~J^Q was usually wider 

than the one for A'^^P whenever a gap can be computed for both. In addition, the gap 
width is often less important than the very ability to compute a gap in the first place. 

For A; = 12 we find similar results for P and Q in Table El with a large number of cases 
in which a gap could be computed for Q, but not for P. The numbers in Table [7] are more 
favorable towards B'^^^Q for the number of detected gaps, whereas for the gap width they 

are more favorable towards A'^^P, although the difference diminishes with diminishing 
ll^iall. For this example we found that the results for Q are not much worse than those 
for BzLQ, which is more costly to compute. 



V 


TOTAL 


TOTAL 


ONLY 


ONLY 


V 


A^'P 


PkpQ 


% 

Gap(B-;2Q) i Gap(^-^P) 


0.25 


40 


69 


2 


31 


0.25 


23 (13) 


19 (9) 


64 


0.125 


42 


161 





119 


0.125 


24 (11) 


22 (10) 


93 





53 


192 





139 





23 (120 


25 (11) 


99 



Gap frequency 



Gap ratio 



Table 5: Gap frequencies and gap ratios for Af^ ^P and B^J^Q when n = 14, m = 25, and 
k = S. 
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V 


P 

TOTAL 


Q 

TOTAL 


P 
ONLY 


Q 
ONLY 


V 


P 


Q 


% 

Gap(Q) i Gap(P) 


0.5 


338 


515 





177 


0.5 


21 (9) 


13 (5) 


15 


0.25 


378 


763 





385 


0.25 


22 (9) 


20 (7) 


56 


0.125 


357 


758 





401 


0.125 


23 (9) 


22 (8) 


74 





366 


757 





391 





23 (9) 


23 (9) 


85 



Gap frequency 



Gap ratio 



Table 6: Gap frequencies and gap ratios for P and Q when n = 14, m = 25, and k = 12. 



V 


TOTAL 


TOTAL 


Al^P 
ONLY 


ONLY 


n 






% 

Gap(B-/2g) i GMAl^P) 


0.5 


759 


766 


11 


18 


0.5 


29 (10) 


15 (5) 


8 


0.25 


794 


880 


1 


87 


0.25 


29 (10) 


21 (7) 


8 


0.125 


767 


859 





92 


0.125 


29 (10) 


24 (7) 


11 





751 


863 





112 





29 (9) 


24 (9) 


15 



Gap frequency 



Gap ratio 



Table 7: Gap frequencies and gap ratios for A^^ and B^j^Q when n = 14, m = 25, and 
k = 12. 



Example 3 

This example is analogous to Example 2. Here, n = 4, ?n = 50, and k = 2. The elements 
of Aj, for j = 0, 1,2,3, have real and imaginary parts that are uniformly distributed in 
[—0.5,0.5], multiplied by a random number which is uniformly distributed in [0,10] for 
j = 0, 1, in [0, 80] for j = 2, and in [0, 77] for j = 3, where rj = 1, 0.5, 0.25, 0. One thousand 
such polynomials were generated for each value of rj. The results are collected in Table [8] 
and Table [9l Again, only the the 1-Norm was used. 

The results are similar to those obtained for Example 2, although here they are signif- 
icantly better for A^^P and B~^J^Q than for P and Q, respectively, which seems to be due 

to the larger value of m. We note that, contrary to the quartic polynomials P and A^^^P, 
Q and B^^^Q are quadratic matrix polynomials for which the roots of fk in Theorem 13.31 
are trivial to compute. 
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V 


P 


Q 


P 


Q 


V 






% 


TOTAL 


TOTAL 


ONLY 


ONLY 


P 


Q 


Gap(Q) i Gap(P) 


1 


22 


25 


8 


11 


1 


15 (7) 


15 (8) 


36 


0.5 


36 


47 


8 


19 


0.5 


14 (7) 


14 (7) 


47 


0.25 


56 


86 


6 


36 


0.25 


14 (7) 


15 (7) 


68 





92 


185 





93 





18 (8) 


17 (8) 


97 



Gap frequency 



Gap ratio 



Table 8: Gap frequencies and gap ratios for P and Q when n = 4, m = 50, and k = 2. 



V 


TOTAL 


PkJiQ 

TOTAL 


ONLY 


Pkl-iQ 
ONLY 


n 


Al'P 




% 

GapCB-i^Q) I Gap(yl-ip) 


1 


272 


124 


162 


14 


1 


18 (8) 


14 (8) 


14 


0.5 


304 


235 


91 


22 


0.5 


20 (8) 


15 (7) 


12 


0.25 


288 


297 


22 


31 


0.25 


20 (8) 


17 (7) 


21 





319 


420 





101 





21 (8) 


20 (8) 


84 



Gap frequency 



Gap ratio 



Table 9: Gap frequencies and gap ratios for A^^^^P and B^J^Q when n = 4, m = 50, and 



k = 2. 



Example 4 

This example illustrates how Lemma 14.11 applied to scalar polynomials (m = 1), can 
sometimes be used to provide gaps between the moduli of zeros of a polynomial when the 
classical scalar version of Pellet's theorem cannot. To this end, we have generated 1000 
scalar polynomials of the form 

p{z) = 2;^° + aigz^^ H h ao , 

where the coefficients Oj are random real numbers, uniformly distributed on [—2, — 1]U[1, 2] 
for j = 3, 5, 11, 13, on [-10, -8] U [8, 10] for j = 4, on [-16, -14] U [14, 16] for j = 12, and 
on [—1, 1] for all J 7^ 3, 4, 5, 11, 12, 13. Theorem 13.31 was applied for = 4, 12. Because the 
matrices involved are all 2 x 2, we have used the 2-Norm, which can easily be computed in 
this case and we have applied Theorem [3]3] only to BJ^^^Q since the matrix multiplications 
are not costly and the results are better than for Q. 

The results are displayed in Table [TU] in the same format as for the previous example. 
The designation A^^P is replaced by p since for the scalar version of Pellet's theorem 
there is no diff^erence between its application to p and to a^^p. The numbers of times that 
a gap could be computed for both = 4 and k = 12 for the same polynomial are listed in 
Table EH 

For A; = 4, more gaps are detected by the classical Pellet theorem than when combining 
Lemma l4. II with the generalized Pellet theorem (Theorem I3.3p . although there are still 21 
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cases where the latter does provide a gap when the classical Pellet theorem cannot and 
that gap is also wider in 29% of the cases when both provide a gap. For k = 12, the 
situation is reversed: using the generalized Pellet theorem for B'^J^Q delivers significantly 
better results than using the scalar Pellet theorem for p. 

As Table [TT] shows. B^^^Q delivered a gap for both k = 4 and /c = 12 in 30 cases more 
than the 4 cases in which the classical Pellet theorem did. 



k 


P 

TOTAL 


TOTAL 


P B^I^Q 
ONLY ONLY 


k 


P 




% 

Gap(B-;2Q) I Gap(p) 


4 
12 


378 
59 


223 
98 


176 21 
9 48 


4 
12 


44 (18) 
26 (13) 


37 (15) 
28 (13) 


29 
82 



Gap frequency 



Gap ratio 



Table 10: Gap frequencies and gap ratios for p and B^J^Q when n = 20, m = 1, and 
k = 4,12. 



p 

TOTAL 



TOTAL 



P 

ONLY 



ONLY 



23 49 4 30 

Table 11: Gap frequencies when n = 20 and m = 1 for both k = A and k = 12. 



Example 5 

For this example, which has the same purpose as Example 4, we have generated two 
classes of scalar lacunary polynomials with random complex coefficients. The first class 
consists of polynomials of the form 

20 7 fi ^ 3 2 

z + ajz + a^z + a^z + a^z + a2Z + aiz + oq , 

where the real and imaginary parts of aj are uniformly distributed on [—2, —1] U [1, 2] for 
j = 1,3,5,7, on [-4,4] for j = 0, on [-10, -8] U [8, 10] for j = 2, and on [-16, -14]U [14, 16] 
for j = 6. 

The second class consists of random polynomials of the form 

+ arz'^ + a^z^ + a^z^ + 032;^ + 02^^ + aiz + oq , 

where the real and imaginary parts of aj are distributed as in the first class. 

The results are presented in Table [12] analogously as for the previous example. Al- 
though the average gap is slightly narrower for BVhQ than for p when k = 2, the reverse 
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is true for k = 6 and there is a clear advantage of B^^^Q over p in the number of computed 
gaps. In addition, there are significantly more cases where using B^^^Q produced a gap 



for both k = 2 and k 
increasingly advantageous. 



6, than when using p. As n increases, using B^J^Q becomes 



n/k 


P 

TOTAL 


TOTAL 


P 

ONLY 


ONLY 


n/k 


P 




% 

Gap(B-;2Q) i Gap(p) 


20/2 


127 


147 


12 


32 


20/2 


53 (20) 


41 (16) 


30 


20/6 


386 


587 


27 


228 


20/6 


31 (13) 


36 (14) 


87 


40/2 


117 


135 


18 


36 


40/2 


49 (20) 


42 (15) 


35 


40/6 


5 


43 





38 


40/6 


19 (16) 


23 (12) 


80 



Gap frequency 



Gap ratio 



Table 12: Gap frequencies and gap ratios for p and B^J^Q when n = 20,40, m = 1, and 
k = 2,6. 



n 


P 

TOTAL 


TOTAL 


P 

ONLY 


ONLY 


20 
40 


97 
3 


123 
24 


14 



40 
21 



Table 13: Gap frequencies when n = 20, 40 and m = 1 for both k = 2 and /c = 6. 



4.2 Scalar polynomials as matrix polynomials 

In the remainder of this work, we concentrate on scalar polynomials. More specifically, 
we propose another way (different from Lemma [4.ip to try and overcome situations where 
the real polynomial fk in Theorem 11.11 does not have positive roots. It is obtained by 
expressing a scalar polynomial as the determinant of a matrix polynomial for which the 
generalized matrix version of the theorem will provide a gap where the scalar version 
could not. As before, it may appear counterintuitive to use the generalized form of the 
theorem to improve the original theorem. There are infinitely many ways to write a scalar 
polynomial in such a form and the best way to proceed is, a priori, not clear. As an 
example, consider the polynomial — + ?>z — 2. It can be written as 

z"^ -y/2 z \ ,/A0\2,/^0l\ , f-y/2 



3 z2 + ^j=d^HU ij^ +il oj^+V-3 V2 
but it can also be written as, e.g., 

3/-2 /-l)=^^*(C ?)"'+(3 o)"+(-2 ~-l 
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Furthermore, one is not limited to using 2x2 matrices, and matrices of larger size can 
also be used. In addition, this approach is, as before, not applicable for every value of k 
in Theorem ll.il On the other hand, it does provide, in general, new upper and/or lower 
bounds on the moduli of the zeros. 

Although it is difficult to predict which of the equivalent matrix polynomial represen- 
tations of a general scalar polynomial will yield the best results, there seems to be some 
indication that better results are obtained if the coefficients are, in a sense, equally dis- 
tributed among the elements of the matrix polynomial. This distribution therefore needs 
to be tailored to the particular polynomial under consideration, which makes it difficult 
to formulate a general method. The following lemma suggests how to proceed by showing 
how to accomplish this with 2x2 matrices for a class of lacunary polynomials. 



Lemma 4.2 Let 



p{z) = az'^ + bz 



where aa ^ 0, and let 



A 



( b 



B 



2^ 



+ 1 



+ az^ -I- /3z -I- 7 , 
4ay 



(7) 



1/2 







c 



V 













D 



2^ 



and W 



+ i{c- — 



1/2 



2^ 



E 



62\ 
4a j 



1/2 







/5 



1/2 














b 






2^/E 


i 


4a 


/3 


-i(^"f - 




l/2\ 


2y/^ 


4a y 


I 



1/2 



^ 4a/ 



/a 

\2^ 

Then the zeros of p are the eigenvalues of the matrix polynomial 

Qeveniz) = ^^"^^ + Bz"/^-^ +Vz + W 

when its degree n is even, and they are the finite eigenvalues of 



(8) 



Qodd{z) = Cz("+^)/2 ^ Dz^n-l)/2 ^ ^^{n-3)/2 ^y^^^r 
when its degree n is odd. 

Proof. The polynomial p can be written as follows: 

p{z) = az"^ + bz""'^ + cz^'"^ + az'^ + (3z + j 

2 



az 



z' + 



.n-2 



a \ z + 



2a 




+ 7 
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Now assume that n is even. After factoring the second order expressions, we obtain 



p{z) 





which is easily seen to be equivalent to det{Qeven{z)), with 

Qeveniz) = Az""^^ + Sz^/^"! +Vz + W , 

where the 2x2 matrices A, B, V, and W are defined in the statement of the lemma. 
When n is odd, we can similarly write p as 

p(z) = det ^ n 1 /9 / 

r- " , 7 - — ^{"-3)/2 _^ -lie 

which is equivalent to det{Qoddiz)), with 

Qodd{z) = Cz("+^)/2 ^ /)^(n-l)/2 + ^^("-3)/2 + + ^ 

where the 2x2 matrices C, D, E, V, and W are defined in the statement of the lemma. 

O 

Lemma 14.21 can now be used in conjunction with Theorem 13.31 and Theorem 13.41 to 
obtain additional and/or improved results for the scalar versions of these theorems, as 
long as the appropriate matrices are nonsingular. In this regard we note that the matrices 
A and V are always nonsingular since we assumed that aa ^ 0, and the matrices C and 
E are always singular. As with Lemma |4. 11 not all of the scalar bounds can be improved, 
but enough can to make it worthwile. When applying a similar process as in Lemma 14.21 
to a general polynomial, the bounds, for a particular value of k, require 0{n) operations 
since the coefficient matrices of the equivalent matrix polynomial are always of size 2x2. 
Similarly to what we had before, the degrees of the new matrix polynoniicils Qeven 

and 

Qodd are half and roughly half that of P, respectively. 

Remarks. Although we will not pursue the matter further, Lemma 14.21 allows for many 
variations on its theme. For instance, a polynomial of the form az'^ + bz"'~^ + cz"'~'^+az+P 
with aa ^ becomes of the form of a polynomial as in Lemma 14.21 with 7 = after 
multiplication by z, so that the lemma becomes applicable to such polynomials as well. 
Moreover, the lacunary polynomial p in Lemma 14.21 is a special case for r = s = 1 of the 
polynomial 

where r > 1 and s > 1 are integers. Lemma 14.21 can easily be modified to include these 
more general polynomials. The polynomials can be further generalized by replacing z with 
a polynomial in z. 
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The following example illustrates how Lemma 14.21 can be used to improve upper and 
lower bounds on the moduli of the zeros of a polynomial and also to detect gaps between 
groups of zeros where the classical scalar Pellet theorem cannot, or to improve the width 
of the gaps that it does detect. 

Example 6 

For this example, we have generated 1000 polynomials of the form ([Tj), with random pa- 
rameters a, b, c, a, 13, and 7 uniformly distributed in [—50,50] and for n = 10,20,40,80. 
Both the upper and lower bounds from Theorem 13.41 and its scalar version were computed, 
as well as the bounds from Theorem 13.31 and its scalar version for k = 2 and k = n — 2, 
i.e., bounds on the gap between 2 and n — 2 zeros. These bounds were computed for both 
p and Qeven- The matrices B and W defining Qeven were nonsingular in all cases. Again, 
as before, because all matrices involved are 2 x 2, we used only the 2-Norm. 

We have computed means and standard deviations of ratios and expressed them as 
percentages and we have listed gap frequencies with the same conventions as in the previous 
examples. The results were collected in the tables below. 

Table [14] lists the ratios of the upper and lower bounds to the moduli of the largest 
and smallest zeros, respectively. The designation "scalar" and "matrix" refers to the use 
of p from ([7]) and Qeven from ([8]), respectively. The last three columns list the percentage 
of cases for which using Qeven yielded a better result than using p for the upper, lower 
and both upper and lower bounds on the moduli of the zeros of p. The upper and lower 
bound ratios using p are quite unaffected by the degree of p, whereas using Qeven yields 
better ratios and smaller standard deviations that improve with increasing degree n. 

Table [15] lists the number of times a gap could be computed with both the scalar and 
generalized (matrix) versions of Theorem l3.3l for k = 2 and k = n — 2. There is a significant 
advantage to using Qeven which, here too, becomes more significant with increasing degree 
n. The same is true for the number of times both gaps could be computed. To pick but 
one instance, for n = 40 and k = n — 2, the classical scalar version of Pellet's theorem 
finds bounds on the gap between the n — 2 smallest and 2 largest zeros 102 times. Using 
the generalized Pellet theorem with Qeven adds another 73 times to that, a significant 
qualitative difference. 

Table [17] contains the means and standard deviations of the gaps and the percentage 
of times Qeven produccd a larger gap than p for those cases in which both produced a gap. 
This was close to half of all cases, while producing a larger average ratio of computed gap 
to exact gap. The gap may be larger for one or the other although the situation could be 
reversed for the bounds defining the computed gap. In other words, Qeven may produce a 
smaller gap than p, but its lower gap bound or its upper gap bound may be better. This 
means that combining the bounds from Qeven and p would improve the gap ratio even 
more. 
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Upper 


Upper 


Lower 


Lower 


% Upper 


% Lower 


% Upper and Lower 


n 


scalar 


matrix 


scalar 


matrix 


better 


better 


better 




10 


119 (26) 


119 (13) 


88 (14) 


83 (11) 


32 


27 


10 




20 


118 (28) 


108 (7) 


90 (15) 


93 (7) 


46 


42 


20 




40 


119 (29) 


103 (4) 


90 (16) 


98 (4) 


65 


59 


39 




QA 

oU 


118 (31) 


1 no /'o\ 
iUz (Z) 


90 (16) 


99 (3) 


72 


71 


51 










Upper and lower bounds for n = 


ID 20 4n RO 








fc = 2 
scalar 




fc = 2 


fc = 2 




fc = n - 2 
scalar 


fc = n - 2 
matrix 


fc = n - 2 


fc = n - 2 


n 


fc = 2 
matrix 


scalar 
only 


matrix 
only 




scalar 
only 


matrix 
only 


10 


41 


25 


21 


5 


47 


26 


26 


5 


20 


81 


73 


30 


22 


66 


61 


31 


26 


40 


123 


161 


19 


57 


102 


157 


18 


73 


80 


127 


192 


15 


80 


147 


221 


11 


85 



Table 15: Gap frequency for n = 10, 20, 40, 80 and k = 2,n-2. 



n 


k — 2 and k — n — 2 
scalar 


fc = 2 and fc = ?i — 2 
matrix 


fc = 2 and fc = ?i — 2 
scalar only 


fc = 2 and fc = ?i — 2 
matrix only 


10 














20 


8 


7 


6 


5 


40 


14 


32 


4 


22 


80 


26 


45 


4 


23 



Table 16: Gap frequency for n = 10, 20, 40, 80 and both k = 2 and k = n — 2. 



n 


fc = 2 


fc = 2 


fc = 2 


fc = n - 2 


fc = n - 2 


fc = n - 2 


scalar 


matrix 


% matrix better 


scalar 


matrix 


% matrix better 


10 


73 (20) 


52 (18) 


10 


70 (25) 


68 (21) 


58 


20 


70 (23) 


66 (19) 


44 


68 (28) 


71 (19) 


46 


40 


74 (26) 


75 (21) 


47 


72 (29) 


80 (19) 


48 


80 


77 (29) 


85 (15) 


46 


71 (31) 


85 (17) 


54 



Table 17: Gap ratio for n = 10, 20, 40, 80 and /c = 2, n - 2. 



Finally, we summarize a few more numerical comparisons, but without a detailed 
breakdown of the results to avoid table overload. 

Applying Lemma |4. 1 1 instead of Lemma 14.21 to the polynomials in Example 6 produces 
results that are generally worse. For k = 2, there are only be about half as many gaps, 
whereas for k = n — 2, the scalar and matrix cases are approximately tied. 

For the odd case of Lemma [4. 21 we have run similar experiments with n = 11, 21, 41, 81 
and k = 2,n — 1 and found that the results for k = 2 are similar, but for k = n — 1, there 
are only about half the number of gaps compared to the even case with k = n — 2. 
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We have also compared the scalar and matrix cases for polynomials of the form az + 
62;^"^ + cz^~^ + az + /3, mentioned in the remarks following Lemma [4.2l for i = 9, 19, 39, 79. 
After multiplication by z, these polynomials are of the form d?]) with n = 10, 20, 40, 80. 
We found that the upper bound and the number of gaps when k = n — 2 were similarly 
improved as for the polynomials in Example 6. However, there was no improvement for 
k = 2. 

Conclusion. We have presented two ideas to overcome some situations where the gener- 
alized Pellet theorem cannot be applied, which at the same time can be used to improve 
bounds resulting from a limit case of this theorem due to Cauchy. Numerical examples 
demonstrate the validity of these ideas. 
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